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Abstract 

A new effect is presented, which changes the emittance during colliding-beam operation in cir- 
cular colliders. If the initial transverse distribution is Gaussian, the collision probability is much 
higher for particles in the core of the beam than in the tails. When small-amplitude particles are 
removed, the remaining ones therefore have a larger transverse emittance. This effect, called core 
depletion, may cause a decrease in luminosity. An approximate analytic model is developed to 
study the effect and benchmarked against a multiparticle tracking simulation. Finally, the time 
evolution of the intensity and emittances of a 208 Pb 82+ bunch in the Large Hadron Collider (LHC) 
at CERN is calculated, taking into account also other processes than collisions. The results show 
that integrated luminosity drops by 3-4% if core depletion is taken into account. It is also found 
that core depletion causes the transverse emittance to be larger when more experiments are active. 
This observation could be checked against experimental data once the LHC is operational. 
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I. INTRODUCTION 



During operation of a circular collider, such as the Large Hadron Collider (LHC) at 
CERN l| , particles are continuously removed or redistributed within the beams by a number 
of different processes, e.g. the collisions, intrabeam scattering (IBS), radiation damping and 
scattering on rest gas. The coupled effect of these processes determine the time evolution of 
the bunch intensities, emittances and the luminosity during a store. 

In this text, the effect of the collisions on the beam distribution is studied using numerical 
parameters for 208 Pb 82+ operation in the LHC. These parameters are given in Table [B If the 
initial transverse bunch distribution is Gaussian, the interaction probability is higher in 
the centre of the bunch. Therefore, the ratio of the number of particles removed in a bunch 
crossing to the initial number of particles over some small transverse distance is much higher 
in the central part of the bunch than in the tails. This leads to a depletion of the core of 
the beam. The transverse emittance of the remaining particles is therefore increasing in the 
absence of damping, which in turn leads to a decreasing luminosity. In principle a similar 
effect is present in the longitudinal plane, since the hourglass effect decreases the collision 
probability for particles ahead of or behind the synchronous particle. This effect is however 
extremely small and can be safely neglected under most realistic machine conditions. 

To study the transverse core depletion effect in detail, an analytic model is developed 
in Sec. [TT] under the assumption that the distribution remains approximately Gaussian but 
changes in size over time. To test the validity of this and other assumptions in the model, it 
is compared with a multiparticle tracking code is described in Sec. IIHI This code includes 
only collisions, synchrotron motion and betatron motion, neglecting other effects such as 
IBS. This somewhat artificial situation is analyzed in order to see the isolated effect of the 
core depletion. 

The analytic model is compared with the tracking in Sec. [IV] and a very good agreement 
is found using numerical parameters for 208 Pb 82+ operation in the LHC. In Sec.|Vl finally, the 
analytic model is extended to include other processes changing the beam distribution. The 
coupled behaviour results in a system of ordinary differential equations (ODEs), which has 
to be solved numerically in analogy with Ref . [2} . Solutions with and without core depletion 
are presented. 
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TABLE I: Nominal parameters for i08 Pb s + ion operation in the LHC, given for the beginning of 
store taken from Ref. [l]. The luminosity reduction factors have been calculated using Eq. (|A4p . 
Three interaction points are considered for ion collisions: Alice (IP2), atlas (IP1) and CMS 

(IPS). 



Parameter 


Value 


Ion species 


208pk82+ 


Beam energy 


2759 GeV/nucleon 


Lorentz factor j re \ 


2963.5 


Bunch intensity iV& 


7 x 10 7 


Bunches per beam 


592 


Normalized transverse rms emittance 


1.5 /im 


Long, emittance at 4 a 


2.5 eV s/charge 


rms bunch length 


7.94 cm 


rms energy spread 


1.1 x 10" 4 


N.o. active interaction points (IPs) 


1-3 


Total interaction cross section a 


515 b 


Optical function f3* y at IP2 


0.5 m 


Optical function /?* at IP1 and IP5 


0.5 m 


Crossing angle at IP2 


70 /irad 


Crossing angle at IP1 and IP5 


285 /xrad 


Geometric luminosity reduction R, IP2 


0.974 


Geometric luminosity reduction R, IP1 and IP5 


0.825 


Peak luminosity 


10 27 cm-V 1 


RF harmonic number h 


35640 


RF gap voltage 


16 MV 
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II. EMITTANCE INCREASE FROM COLLISIONS 



To derive an approximate expression for the increase in emittance caused by core deple- 
tion, the distributions of the betatron action of the incoming and outcoming bunches in a 
collision are calculated. To highlight the features of the process, first the simplified case with 
no other processes acting on the beam is studied analytically. The formulas are derived for 
unequal beams, and then simplified to the case of equal beams. Numerical models including 
other processes, such as IBS and radiation damping, are discussed in Sec. IVl 

The model is constructed under the approximation that the transverse distributions re- 
main close to Gaussian during the whole store. Furthermore, it is assumed that the hourglass 
effect and crossing angle have a negligible influence on the shape of the distribution of the 
colliding particles, which makes the integrals analytically solvable. The angle and hourglass 
effect are instead included in an approximate way once the distribution is known. These 
approximations are justified in Sec. IIV1 

The total number of particles removed, C sc , per interaction cross section a during a single 
bunch crossing is given by an overlap integral of the densities of the two bunches 3|, |4j: 

C SC = MN 1 N 2 J p 1 (x,y,s,T)p 2 (x,y,s,T)dxdydsdT, (1) 

where M = (vi — v 2 ) 2 — (vi x v 2 ) 2 / c 2 is a kinematic factor . The two bunches are 
assumed to move with opposite velocities Vi so their centres have the longitudinal coordinate 
s = ±vr. Both centres are at s = at the interaction point (IP) at time r = 0. Furthermore, 
Pi(x,y,Zi) is the density of bunch % normalized to one, Ni its intensity, and (x, y) are the 
transverse coordinates. All integrations are to be carried out on the interval [— oo, +oo] and 
this convention holds for all subsequent integrals unless indicated otherwise. 

To calculate the distribution in betatron action, Eq. ([1]) has to be generalized to include 
the angular distributions of the bunches. It is assumed that the distributions in the three 
planes are independent and can be decoupled, which means that the crossing angle = 
and the hourglass effect is weak. Then we have V\ = —v 2 and the transverse coordinates 
axes are equal for both beams. With \vi\ = v the kinematic pre-factor becomes M = 2v and 
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C xr can be written as 



C sc = 2vN l N 2 / pi x {x, x[)pi y (y, y[)piz{s - vr) x p 2x (x, x' 2 )p 2y {y, y 2 )p2z{s + vr) 

dx dy dx[ dx' 2 dy[ dy' 2 ds dr. (2) 



3* y ( u 2 + (a* xy u + f3* xy u') 2 



The incoming bunches are assumed to be Gaussian: 

3* 

p iu {u,u ) = exp 

Zna ui 

. 1 
p iz (s±vT) = —^= — exp 



2^ 
s ± vt) 2 



(3) 



for u = x,y and i = 1, 2 for the two beams. Here a ui = a/ '3* y e U i are the transverse beam 
sizes, e U i are the transverse emittances, <7i Z is the rms bunch length, and (3* , a* y ) are the 
optical parameters at the IP, which are assumed to be equal for both beams and in both 
planes. 

Using Eq. ([3]), all integrations in Eq. ([2]) are now carried out, except over x,x[. The 
remaining integrand A gives the number of reactions per cross section for particles in bunch 1 
in a phase space element dxdx[ during a single bunch crossing. Completely analogous 
calculations can be carried out in the vertical plane. 

Normalizing by the total number of interactions, the distribution of collision points is 
thus 

2xx 1 a xy 3 xy + x 1 8 xy + x 



x PxyV °x\ + °x2 

X = —z- — exp 

Z7ra xl a x2 



1 + a* 2 + ^ 



2a 2 xl 




(4) 



For equal beams (a ul = a u2 ) and a* = 0, A simplifies to the result in Ref. [6j, where it is 



shown that the transverse distribution of the collision points in that case is narrower than 
the incoming bunch by a factor 

Changing to action-angle variables ( J x , <p x ) through 



'2J x 3* y cos 



^ (sin <f) x + a* xy cos <p x ) 



(5) 



and averaging over <p x , results in the distribution Aj of the betatron action of the colliding 
particles: 

— 2— /i r 9 i 



Here Iq is a modified Bessel function. 

Assuming that m particles are removed in total from the bunch, the number of par- 
ticles left in a small element dJ x after the crossing is (Nipj — m\j)dJ x , where pj = 
exp (— J x /e x \) /e x i is the incoming distribution of J x obtained by integrating Eq. Q over 
(j) x . The density p rem of the remaining particles, normalized to unity, is thus 

PremK-'x) — jy _ m ' \ ' > 

The expectation value of J x of the incoming bunch is e x \, while for the outcoming bunch 
after the collision it is 

poo 

tlx — / Jx Prem (J x )dJ x = 6,1(1 + (8) 

JO 



where 



2 



maz 



C 2{N l -m){al 1 + al 2 y (9) 
Thus, the change in emittance during the crossing is de x i = ei x — e^i = ( e xi- Averaged over 
one turn, the emittance blowup per time becomes 

= e ^lC/rev^IP = TjT, (10) 

where / rev is the revolution frequency and nrp the number of IPs. 

The number of removed particles is given by m = C sc a, inserting Eq. ([3]) in Eq. ([2]) 
and carrying out all integrations. So far, the luminosity reduction factor R, including the 
hourglass effect and the crossing angle, has been neglected. An approximate way of including 
it, which according to comparison with the tracking simulation in Sec. [IV] is accurate, is to 
use the distribution Aj calculated above but include R in the calculation of m. We therefore 
have 

iViiV 2 C7 . , 

m = aR . (11) 

27ryV*i + ^)(^i + ^ 2 ) 

A general expression for R is given in Appendix [Aj 

With Eq. pip , the rise time T c can be expressed in known parameters. Normally m<iVj, 

so for simplicity only first order in m/iVj is kept. This gives 



4vr + al 2 f 2 Jo* + al 



T c = ¥ - — (12) 

N 2 o-Ra 2 xl f rew n IP 
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An analogous expression holds for the rise time of the vertical emittance. If the beams are 
round, the rise time of beam i can be written as 

_ 4^/%, (e xyi + e xy ,f 2 

c ~~ / AT D~f ' V / 

iv janj TCV nip 

where the beam sizes have been written in terms of the transverse emittances e xy i. Here 
i = 1, 2 and j = 2 for i = 1 and vice versa for the two beams. If the beams are equal, we 
obtain finally 

c NaRf iev n 1P ' 1 1 

This is the average time rate of the emittance change over the first turn for an initially 

Gaussian bunch and Eq. (j7j) gives the exact non-Gaussian distribution of the outcoming 

bunch (neglecting crossing angle and hourglass effect). This distribution then turns in phase 

space during one revolution before it enters the IP again, and to obtain an exact distribution 

at later turns the integrations leading to Eqs. (jlj)-((7j) should be repeated recursively. 

The tracking shows, however, that for the LHC parameters the perturbation from a 
Gaussian is very small (see Sec. IIVj) . Therefore, Eq. ffl4|) can be used also at later times 
without significant loss in accuracy. Another necessary condition for this approximation is 
that m < iVi, so that the bunch remains matched after the crossing. Otherwise the rotation 
in the transverse phase space has to be taken into account. 

A simple model of the coupled time evolution of the bunch intensity and emittance can 



now be constructed. With equal beams, the instantaneous luminosity C is j^] 

£Mt) 3fe< < i5) 

where is the number of bunches. The rate of removal of particles at an IP is given by 

aC, so with Eq. (1151) we have (for one bunch colliding at nip IPs): 

dN(t) = N 2 (t)f iev n JP _m_ 

dt 1 1 47re xy (t)f3* y T c (t) { ' 

Here Tc is defined as the instantaneous lifetime due to collisions. The time evolution is given 

by Eq. (fl6l) coupled with Eq. (flQl) . which in the case of equal beams simplifies to 

de xy = I6ire xy (t)p* y 

dt ~ N(t)aR(t)f iev m P ' ^ ' 

It should be kept in mind that, when a strong betatron coupling is present, the emittance 
increase is shared between the planes but because of the analogous emittance increase in 
the vertical plane, the net effect for each plane is still given by Eq. (TfJ 



In the general case, when a crossing angle is present, R is a function of the emittances 
and therefore time dependent. The system of ODEs then has to be solved numerically. If 
there is no crossing angle, or when its effect is very small, R is constant and the ODEs have 
the analytic solution 

N(0) 



N(t) 



1 + 4*- 



4/5 



e xy {t) = e xy (0) (l + ^-J 



1/5 



(18) 



where T c0 = T c at t = 0. This solution can be compared to the case when core depletion is 
not taken into account. Then Eq. (fT6l) can be solved directly, assuming a constant emittance, 
to yield 



N(t) 



N(0) 
1 + f- 

IcO 



Using Eq. ffTo^) . the luminosity at one IP with core depletion is 

f mv k h N 2 (0)R 4A; 6 iV(0) 1 



and 



C(t) 



C(t) 



Wo)^(i + %) Tconipa 

f mv k b N 2 (0)R 4A; fe iV(0) 



1 + 



5f 



9/5 



4^(0)/?*, (l + ^) Tc ° nipa 
without it. Finally, the integrated luminosity including core depletion is given by 



(19) 



(20) 



(21) 



C(t) dr 



fcbiV(O) 
n IP o" 



/ 



V 



5t 



4/5 



(22) 



while if the core depletion is neglected, it reduces to 



C(t) dr 



fcbiV(O) 
nipcx 



(23) 



With the parameters given in Tabled! we have T c0 = 29.9 h for n^> = 3. This is compared 
with the strength of other effects in Sec. [V] 



III. MULTIPARTICLE TRACKING 



The tracking simulation program follows two bunches containing a number of macro 
particles. In the simulations presented here, 5 x 10 4 particles were used to represent a 
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FIG. 1: Schematic illustration of a collision between two bunches at an IP. The directions of the 
movement are indicated by the red arrows in the front of the bunches. The y-coordinates coincide 
(yi = 2/2)- All distances are measured in the lab frame. 



nominal bunch of 7 x 10 7 208 Pb 82+ ions in the LHC. The 6D coordinates of the particles are 
updated on a turn-by-turn basis by three routines: betatron motion, synchrotron motion 
and collisions. All other processes are neglected in order to isolate the core depletion effect. 

Betatron motion is represented by a rotation in normalized phase space by an angle given 
by the machine tune and chromaticity. Synchrotron motion is implemented by a change in 
energy and longitudinal momentum. Both routines are taken from Ref. 7j], where a more 
detailed description can be found. 

To simulate the collisions, the program loops through all particles and calculates for each 
of them an interaction probability P\ as a function of its coordinates and the distribution of 
the opposing bunch. A random number is then sampled to determine if an interaction takes 
place, in which case the particle is removed. To calculate Pi, the movement of a particle in 
bunch 1 through bunch 2 at an IP is considered. If bunch 1 contains only one particle, we 
have Pi = crC sc , with C sc given by Eq. (CD) without approximations. 

In the general case, a crossing angle <fi 7^ has to be taken into account, which is assumed 
to be in the horizontal plane. The different coordinate systems used are defined in Fig. [T], 
where the s$ axes are fixed and the Zj axes move with each bunch. For simplicity, both 
densities will be expressed in the S2-X2 system. Since a negligible transverse magnetic field 
is assumed at the IP (the experimental chambers are usually constructed in such a way that 
this is fulfilled) a particle in bunch 1 with spatial coordinates (xio, yw, z w ) at time r = 
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and transverse angles (x[,y[) follows approximately a straight line given by 



Xi 






— X 1q Z\q 




x 10 


yi 




2/io 


- 2/io^io 


+ si 


y'w 


S\ 











1 



(24) 



where it is assumed that the longitudinal coordinate changes in time as 

si = z 10 + vr. 

In the rotated system of bunch 2, Eq. (|24p is transformed to 



(25) 



x 2 




c 





-s 




Xi 




^20 




X 20 


2/2 







1 







2/i 




2/20 


+ S 2 


2/20 


$2 




S 





c 




Sl 









1 



(26) 



where C = cos 4> and S = sin <p. Eq. ff26l) contains five unknowns and three equations. The 
last equation can be used to express Si in terms of s^, which upon insertion in the first two 
equations can be used to identify the remaining coefficients. The solution is 



2^20 — (x\o — x' w Ziq) / A 

X 20 = { X 10^ ~ 

2/20 = yio-y[ (z 10 c + x 10 S)/A 

2/20 = y'w/ A 

s 2 (t) = (vr + z w )C + (xio + vrx' w )S (27) 

where A = C + x' 1Q S and Eq. ff25l) has been used to derive the r-dependence of s 2 . 

The density function p\ for a single particle, needed in Eq. flTJ, can be modelled by the 
Dirac 5-function: 

Pi(x2,y2,s 2 ,r) = 5(x 2 - [x 20 + x' 20 s 2 })5(y2 - [y 2 o + 2/ 2 o s ])^( s 2 - s 2 (r)). (28) 

The particles in the opposing bunch are sorted in discrete bins along the directions (x, y, z 2 ) 
in order to obtain the density p 2 . It is assumed that the transverse distributions are inde- 
pendent around the IP, and that the longitudinal density p 2z does not depend on x or y. The 
transverse binnings are performed using a constant (3 xy = (3* , but it has to be accounted 
for that the distribution of bunch 2 changes along s 2 with f3 xy (s2), given by 



Acta) = PL i + 



(29) 



XI) 
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This can be modelled through a parameter k(s2), which gives the ratio of the width of 
bunch 2 at S2 7^ to the width at $2 = 0: 

The widening of the beam is thus expressed as P2 X (%2, $2) = P2x(x 2 /k, 0)/k = p2 X (x 2 / k)/ k. 

The interaction probability Pi = aC sc for a particle is then obtained by integrating 
Eq. (P) with pi given by Eq. (1281) . All integrations except over S2 can be carried out directly. 
Using that the kinematic pre-factor simplifies to M = 2v cos 2 (0/2), the result is 

_ 2iWxcos 2 f f PU^^Pjyi^^) /r 11 Z W C + X 10 S \ ■ 

(31) 

Here pl x (x), p2 y {y) are the transverse densities of bunch 2 at the IP. The integral in Eq. fl3Tl) . 
which is solved on every turn for every particle, is replaced in the code by a sum over all 
bins that the particle passes through. 

Using this mathematical model for the collisions, the core depletion effect as well as the 
hourglass effect are automatically accounted for. Simulation results from the tracking code 
are presented together with results from the analytic model in Sec. IIVI 

IV. SIMULATION RESULTS WITH COLLISIONS ONLY 

For the comparison between the analytic model in Sec. [Til and the tracking described in 
Sec. Hm two cases are considered: Either the crossing angle is = at all IPs, or <ft is given 
by Table [TJ All other parameters are taken from Table [TJ in both cases and all three IPs are 
assumed active. All other processes except collisions are neglected. 

The results for the first case are shown in Fig. [21 where = implies that R is constant 
so that the analytic ODE solutions are valid. As expected, it can clearly be seen that there 
is an emittance increase in the tracking, which is arising solely from the variation in collision 
probability between the core and the tails of the beam. In this example, when other effects 
are not taken into account, the effective emittance increase is around 20% over 10 h. 

The agreement in bunch population and luminosity between the analytic model with 
core depletion included and the tracking is excellent. The neglection of the core depletion 
introduces a small error, which corresponds to a 5% change in integrated luminosity during 
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10 



t(h) 



t(h) 

FIG. 2: The time evolution of the bunch intensity N, the transverse emittance e xy and the 
luminosity C for 208 Pb 82+ ions in the LHC, with three active IPs, when all processes but the 
collision are neglected, as obtained with the tracking and the analytic solution of the ODEs with 
and without core depletion. The crossing angle is set to (ft = at all IPs and all other parameters 

are given in Tab. HI 



a 10 h store. A good agreement between the methods is also found for the transverse 
emittance, although a small difference can be seen towards the end of the store. 

The analogous results for the second case with a non-zero crossing angle are shown in 
Fig. [3l The tracking is compared both with a numeric integration of the ODEs using 
Mathematica s| (with R given by Eq. ( 1A4I) evaluated at every integration step) and the 
analytic solution with an assumed constant R given in Table [B In this case, the analytic 
solution is very accurate as R does not change significantly. Again, an excellent agreement 
is found in luminosity and bunch intensity, while there is a small discrepancy in emittance. 

A closer examination of the transverse profiles shows that the discrepancy comes from 
the approximation of Gaussian bunches. The distributions in the tracking are not strictly 
Gaussian, but very similar to a Gaussian with a larger standard deviation, which causes 
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FIG. 3: The time evolution of the bunch intensity N, the transverse emittance e xy and the 
luminosity C for 208 Pb 82+ ions at IP2 in the LHC, with three active IPs, when all processes but 
the collision are neglected. Results obtained by the tracking simulation described in Sec. IIIH the 
analytic model including core depletion (analytic and numeric solution), and the analytic model 
neglecting core depletion, given by Eqs. (|19p and (|2ip . The numeric parameters used are given in 
Tab. [U except the crossing angle which has been set to zero. 



a small variation in emittance even though the luminosity and bunch population agree. 
Fig. H] shows the bunch profile from tracking and the analytic Gaussian distribution for the 
case with 0^0. The curves are snapshots of the distribution at 2.5 h intervals, with the 
uppermost curve corresponding to t — h. A small difference can be seen in the tails. 

In spite of the small discrepancy in emittance, the similarity between the two distributions 
throughout the store and the excellent agreement in luminosity and bunch intensity in Figs. [2] 
and [3] justify the approximations in the analytic model. 
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FIG. 4: The transverse bunch file at different times as simulated by tracking and an analytic 
model with a Gaussian with the emittance given by Eq. (|18p . The top curve is at t = and each 
lower curve corresponds to a time 2.5 h later. The distributions are normalized to a /^-function of 
67.5 m and the integral of each curve gives the total bunch population at that time. The crossing 

angles are given in Tab. [H 

V. LUMINOSITY TIME EVOLUTION INCLUDING OTHER PROCESSES 



To model a real machine, other effects such as IBS and radiation damping have to be taken 
into account and the longitudinal emittance en has to be introduced as another dynamic 
variable. In Ref . [2j , the time evolution of the bunch population, emittance and luminosity 
was calculated through numerical solution of a system of coupled ODEs. In this section, an 
analogous calculation is carried out, both with and without the inclusion of core depletion. 

The time evolution of the emittances and the bunch intensity can be described by the 
following system of six ODEs with i — 1, 2 (expanding the result in Q): 



-xyi 

dt 
de h 
dt 
dNi 
~df 



'-xyi '-xyi 



+ 



-xyi 



+ 



-xyi 



MCS 



txyi j 



N; 



i] i ^-xyij ) ^lij 



■ rad,z 

Ni_ 
T 



(32) 



Here the following notation has been introduced: T rax x iXy , T rax x tZ are the radiation damping 
times in the transverse and longitudinal planes, Tj^s, X y, ^ibs,« are the emittance rise times 
due to IBS, Tmcs is the rise time due to multiple Coulomb scattering on rest gas, and T gas 
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is the lifetime caused by inelastic scattering on rest gas. In the LHC, quantum excitation 
is too weak to have an influence and is therefore neglected. To solve Eqs. fl32l) numerically, 
Mathematica was used, taking into account the different crossing angles and /3* y at the IPs 
shown in Table [H 

The IBS rise times were calculated with mad-x, where a generalized version of the 
Bjorken-Mtingwa model was used 0, [h], ll|. The evaluation of Tms^ and T\bs,x V is done 
off-line on a grid of points and interpolated at run-time as in Ref. 2j. Radiation damping 
times, as well as T MCS and T gas , are calculated using standard formulas [4|. In order to show 
the strengths of the different processes, Table |V] presents numerical values of the lifetimes 
and rise times in the beginning of the store using the starting parameters in Table [B As can 
be seen, both processes related to the rest gas are negligible. 



TABLE II: Initial values of the rise times and lifetimes resulting from different processes, 

calculated using the values in Table [H 





Lifetimes 


Rise times and damping times 




TjC ^gas 


T\BS,xy Tibsj T c0 T raA>xy T Tad>z Tmcs 


nip = 1 
nip = 3 


22.4 647.5 
7.5 647.5 


13.2 7.8 89.7 12.6 6.3 43598 
13.2 7.8 29.9 12.6 6.3 43598 



Eqs. ( 1321) are based on the assumption that the beams remain Gaussian throughout the 
whole store. In particular, the expressions for Tc, T c and the IBS rise times are only valid for 
this case. It was shown in Sec. |IV]that the collisions themselves only cause small deviations 



from a Gaussian distribution and measurements at RHIC [12] have shown similar results 
for IBS. Furthermore, radiation damping can be represented by a multiplication of the 
oscillation amplitudes by a decay coefficient, which is the same for all particles. Therefore, 
this does not change the shape of the distribution, only the standard deviation. Finally, the 
beam-gas scattering processes are too weak to have any significant influence. Therefore, it 
is a fair approximation to assume that the beams keep their Gaussian shape throughout the 
store. 

In the case of equal beams, which is studied here, the system f )32|) contain only three 
equations. The last term in the first equation represents the core depletion and to study 
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FIG. 5: The time evolution, given by Eqs. (j32[) . of the bunch intensity, luminosity at IP2, 
transverse rms emittance and longitudinal emittance during a 10 h store at top energy with 
colliding 208 Pb 82+ beams in the LHC. Results are shown from calculations with (dotted lines) and 
without (solid lines) core depletion for the cases of one (green lines), two (red lines) or three 
(black lines) active IPs taking collisions. When only one IP is active, IP2 has been chosen. 



its effect Eqs. (I3"2"j) were solved also with this term excluded for nrp = 1, 2, 3. The result- 
ing luminosity, bunch intensity and emittances are shown in Fig. As can be seen, the 
emittance is shrinking, since radiation damping is stronger than IBS. There is a small but 
notable difference between the luminosity with and without core depletion. The ratio of the 
integrated luminosity over 10 h with core depletion included to the case without it is 0.97 
for one active IP and 0.96 for three. 

A striking difference can be seen in the qualitative behaviour of the time evolution of the 
transverse emittance in Fig. Without core depletion, the emittance shrinks faster with 
more IPs active, since more particles are removed through collisions in this case. Therefore 
the effect of IBS becomes weaker with time while radiation damping is independent of the 
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intensity. When core depletion is included, the emittance shrinks instead faster with only one 
active IP, as T c scales linearly with n^p as shown in Eq. (JHJ). This qualitative behaviour is an 
important observation which could be checked experimentally once the LHC is operational. 

The longitudinal emittance is also shrinking, but in this case the emittance shrinks faster 
when core depletion is taken into account. This can be understood by considering that IBS 
is weaker in the longitudinal plane when the transverse emittances are larger. 

VI. CONCLUSIONS 

A new effect that increases the emittance in circular colliders with Gaussian beam profiles 
has been presented. Since the interaction probability in the collisions is much higher in the 
centre of the bunch than in the tails, the core of the beam is depleted, so that the emittance 
of the surviving particles is larger. 

The effect has been studied first through a simple multiparticle tracking simulation, which 
makes no assumptions on the shape of the beam distribution. The results show that the 
emittance is indeed growing due to the collisions and that the transverse distribution remains 
close to Gaussian. 

To describe the effect analytically, the expectation value of the betatron action was cal- 
culated before and after a bunch crossing with an initially Gaussian bunch. The emittance 
increase was averaged over time to form an ODE, which coupled with another ODE describ- 
ing the evolution of the bunch intensity describes the time evolution under the assumption 
that the distribution remains Gaussian. Results from the analytic model agree very well 
with the tracking. 

Finally a system of ODEs including also other effects changing the beam distribution and 
intensity, such as IBS, radiation damping and scattering on rest gas, was solved numerically 
for the case of 208 Pb 82+ operation in the LHC. It was shown that a qualitatively different 
behaviour of the transverse emittance is expected when core depletion is included: The 
emittance is larger when more IPs are active, as opposed to the expected behaviour without 
core depletion. This prediction could be verified experimentally. Quantitatively, the core 
depletion effect is expected to introduce corrections of 3-4% to the existing calculations of 
integrated luminosity for 208 Pb 82+ beams in the LHC. 
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APPENDIX A: LUMINOSITY REDUCTION FACTOR 



Reductions to the luminosity due to the hourglass effect [13] and a non-zero crossing 



angle [l^J] are well-known phenomena. For the purposes of this text, the results from Ref. 
will be rewritten on a slightly different form. 

To obtain a general reduction factor for the luminosity with Gaussian bunches, relative to 
the limiting case without crossing angle and hourglass effect, we start from the total number 
of interactions in a bunch crossing given by Eq. ([1]). Using the coordinate systems defined 
in Fig. [1], the beam distributions in the Zi-system moving with bunch i can be written as 

y2 \ / ,,2 \ / „2 

» i 



27T(Txi yl'KGyi \/lTXO z 



All integrations are carried out in the fixed x — s system, using the transformation 



Xi = x cos ipi + s sin 

Vi = y 

Si = —x smipi + s cos (A2) 

where ipi = 0/2 and ip2 = —0/2. Inserting the transformed distributions in Eq. ([T]), the t, 
x, and ^/-coordinates can be integrated to yield 

C sc = R, (A3) 

2 ^PxyV ( e xi + e x2 ){e y i + e y2 ) 

where R is a total reduction factor coming from the crossing angle and the hourglass effect, 
which, using Eq. ( |29l) . is given by 



R 



I I 9 2 cos 2 ipi 1 sin 2 i/>i 
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■da. (A4) 



In the case of equal beams, Eq. (1A3I) is equivalent to Eq. (2) in Ref. [14|. The integral 
in Eq. (IA4I) is not analytically solvable in the general case and was therefore integrated 
numerically throughout this text. With zero crossing angle the integral simplifies to the 
results in Ref. [13( and if the hourglass effect is neglected (short bunches) the reduction 
factor with finite crossing angle in Ref. [ijj] is obtained. 
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